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Ph' Abstract 



This work is a continuation of our previous work (JMP, 48, 12, pp. 122103-1-122103- 
20, 2007), where we constructed the non-relativistic Lee model in three dimensional Rie- 
mannian manifolds. Here we renormalize the two dimensional version by using the same 
methods and the results are shortly given since the calculations are basically the same as 
in the three dimensional model. We also show that the ground state energy is bounded 
■ from below due to the upper bound of the heat kernel for compact and Cartan-Hadamard 

manifolds. In contrast to the construction of the model and the proof of the lower bound 
of the ground state energy, the mean field approximation to the two dimensional model is 
not similar to the one in three dimensions and it requires a deeper analysis, which is the 
main result of this paper. 



PACS numbers: ll.10.Gh, 03.65.-w, 03.65.Ge 

1 Introduction 

The Lee model was originally introduced in [1] as an exactly soluble (in principle eigenstates 
and eigenvalues can be exactly found) and a renormalizable model that describes the interaction 
between a relativistic neutral bosonic field "pions" and two neutral fermionic fields "nucleons" . 
It is assumed that the nucleon can exist in two different intrinsic states. The particle corre- 
sponding to the Bose field is called 9 and the particles corresponding to the intrinsic states of 
the nucleon are called V and N particles. The fermionic field corresponding to the V and TV 
particles are assumed to be spinless for simplicity. Only allowable process is given by 

v — N + e (i) 

and the following process is not allowed 

N ±=?V + 6 , (2) 

which makes the model rather simple. Although this model is not realistic, the important 
features of nucleon-pion system can be understood in a relatively simple way and one can get 



I 



rid of the infinities without applying perturbation theory techniques. Moreover, the complete 
non-relativistic version of this model that describes one heavy particle sitting at some fixed point 
interacting with a field of non-relativistic bosons is as important as its relativistic counterpart. 
It is much simpler than its relativistic version because only an additive renormalization of the 
mass difference of the fermions is necessary. It has been studied in a textbook by Henley 
and Thirring for small number of bosons from the point of view of scattering matrix [2] and 
there are various other approaches to the model [31 EU EJ EJ [7J El E]. It is possible to look at 
the same problem from the point of view of the resolvent of the Hamiltonian in a Fock space 
formalism with arbitrary number of bosons (in fact there is a conserved quantity which allows 
us to restrict the problem to the direct sum of n and n + 1 boson sectors). This is achieved 
in a very interesting unpublished paper by S. G. Rajeev [ID], in which a new non-perturbative 
formulation of renormalization for some models with contact interactions has been proposed. 

This paper is a natural continuation of our previous work [XT] and we extend the three 
dimensional model constructed there to the two dimensional one. In [11], we discussed the 
non-relativistic Lee model on three dimensional Riemannian manifolds by following [10] and 
renormalized the model by the help of heat kernel techniques with the hope that one may 
understand the nature of renormalization on general curved spaces better. In fact, the idea 
developed in [TO] has also been applied to point interactions in [121 EU [33] and to the relativistic 
Lee model in [15]. In this paper, we are not going to review the ideas developed in [10] and 
|llj . Instead, we recommend the reader to read through these papers. The construction of 
the model is exactly the same as the one in three dimensions. It is again based on finding 
the resolvent of the regularized Hamiltonian H e and showing that a well-definite finite limit of 
the resolvent exists as e — > + (called renormalization) with the help of heat kernel. We then 
prove that the ground state energy for a fixed number of bosons is bounded from below, using 
the upper bound estimates of heat kernel for some classes of Riemannian manifolds. Finally, 
we study the model in the mean field approximation for compact and non-compact manifolds 
separately. Although the construction of the model and lower bound of the ground state energy 
are based on the same calculations as in three dimensions, the mean field approximation in two 
dimensions requires a deeper analysis as we shall see, which is the main result of the present 
work. 

The paper is organized as follows. In the first part, we give a short construction of the 
model and show that the renormalization can be accomplished on two dimensional Riemannian 
manifolds. Then, we prove that there exists a lower bound on the ground state energy. Finally, 
the model is examined in the mean field approximation. In appendix, we prove an inequality 
which we use in the mean field approximation. 

2 Construction of the Model in Two Dimensions 

In this section, we give a brief summary of the construction of the model in two dimensions since 
it has been basically done in [11] for the three dimensional model. We start with the regularized 
Hamiltonian of the non-relativistic Lee model on two dimensional Riemannian manifold (Ai, g) 
with a cut-off e. In natural units (h = c = 1), the regularized Hamiltonian on the local 
coordinates x = (xi,x 2 ) G M. is 

H* = H + H I>e , (3) 



where 

Ho = J m d 2 g x <j>l(x) i-^l + m) cj> g {x) , (4) 

Hl, e = /i(e) y—J~^J + X J d l X K e{ x , «! flO (09^) ® — + 4>l(x) • ( 5 ) 

Here <f) g {x), <f) g (x) is the bosonic creation-annihilation operators defined on the manifold with 
the metric structure g and A is the coupling constant, and x, y refers to points on the manifold 
j\A. Also, K e (x, a; g) is the heat kernel on a Riemannian manifold with metric structure g and 
it converges to the Dirac delta function S g (x, a) around the point a on Ai as we take the limit 
e — > + . For simplicity, we have changed the notation for the heat kernel to K s (x, y; g) instead 
of writing K s / 2m (x, y; g) which was used in [11]. We also assume stochastic completeness, that 
is 



/ 

JM 



d g x K s (x,y;g) = 1 . (6) 



<M 

Similar to the flat case, /j,(e) is defined as a bare mass difference between the V particle (neutron) 
and the iV particle (proton). Although the number of bosons is not conserved in the model, 
one can derive from the equations of motion that there exists a conserved quantity 



Q = -[ ^ ) + / <x <P\{x)<P g {x) 



■IM 



2. 



Therefore, we can express the regularized Hamiltonian cts di> 2 x 2 block split according to C 

fje_ F= f H -E A J M d 2 g x K e (x,a;g)(f>l(x) \ 

\Xj M d 2 g xK e (x,a;g)cf> g (x) H - E + ^(e) )' { ' ] 

Then, the regularized resolvent of this Hamiltonian in two dimensions is 

1 / a e (3\ \ ( a e b\ N 



where 



*W = H^E = ft Z J = I i ) < 8 > 



b\ Q-\E) b € 



Hg — E Hq — E Hq — E 

5 e = <5>;\E) 

b e = A / d 2 g x K e (x,a;g)(j) g (x) . (9) 
jm 

Here E should be considered as a complex variable. Most importantly, the operator § t (E), 
called principal operator, is given as 

$ e (E) = H - E + Me) - A 2 f d 2 g x d 2 g y K e (x, a;g)K e (y, a; g) 3 (x)— 1— <j>\{y) . (10) 

Jm 2 -no - E 

After performing normal ordering of this operator, we get 

$ e (£) =H -E- A 2 [ ds [ d 2 g xd 2 g y K s (x,a;g)K 3 (y,a;g) 

Je/2 Jm 2 



/oo 
ds K s {a, a; g) [e~< m ^ - e -(s-e)(H +m-E)^ _ (n) 

If we choose fi(e) as 

POO 

/i(e)=/i + A 2 / ds K s (a,a;g)e- S(m -> 1) , (12) 



where /i is defined as the physical energy of the composite state which consists of a boson and 
the attractive heavy neutron at the center a, and take the limit e — > + , we obtain 



/"OO 

H -E + ii + \ 2 ds K s (a, a; #) [e^" 1 "^ - e -^ + m -E)^ 
Jo 



poo p 

-A 2 / ds d 2 g xd 2 g y K s (x,a;g)K s (y,a;g)<l>l(x)e- s ^ +2m - E U g (y). (13) 
Jo Jm 2 

This is the renormalized form of the principal operator so that we have a well-defined explicit 
formula for the resolvent of the Hamiltonian in terms of the inverse of the principal operator 
^~ 1 (E). As we will elaborate in more detail later on, the bound states arise from the roots of 
the equation 

$(£01*) = 0, (14) 

corresponding to the poles in the resolvent. Hence, the principal operator determines the bound 
state spectrum. 



3 A Lower Bound on the Ground State Energy 



Following exactly the same method developed for three dimensions [TT] , we again find the upper 
bound of the norm of the operator U'{E) as 



^ 2 r(2) r ,^,-si nm+ ,- E) r A .. r A .. f , 

(ttiM 2 ) i/2 



\u\E)\\ ^ n f{Y^ 1°° ds s e ~ s{nm+fl ~ E) I d ^J o du *l du 



* [K2s(u 1+ u 3 )(a, a; g)] 1/2 [^s^+u^ia, a; g)] 1/2 , (15) 

where T denotes the gamma function. For each class of manifolds, there are different upper 
bounds on the heat kernel so we will consider them separately. We will first consider Cartan- 
Hadamard manifolds. The diagonal upper bound of the heat kernel for two dimensional Cartan- 
Hadamard manifolds is given as [161 EH] , 

K ^9)<J^y (16) 

for all x G Ai, s > 0, and C is a constant. Using ffTBT) and performing u 3 integral in (|15|) . we 

get 

\2 TVO\ POO 

\\U'(E)\\ < n Cm , \{ / ds e -^m+n-E) 
M I )\\- r(l/2) 2 7 

pi 1 1 Pl-m 1 j 

X k dUl ^(l-«i)^ J dU2 ^(l-n 2 )V 2 
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< n Cm 



A 2 T(2) 



ds e - s i nm +^~ E ) 



du 



1 



1 



-1 2 



MV2(1- U )V2 



r(i/2) 2 

Evaluating the integrals give the following result 

\\U'(E)\\ < n Cmn A 2 T(2)(nm + fx — Ey 1 . 
Then, the inequality above implies a lower bound for the ground state energy 

E gr > nm + u — nC\ 2 m , 

where 



(17) 



(19) 
(20) 



c = cvr(2) . 

The diagonal upper bound of the heat kernel for two dimensional compact manifolds with Ricci 
curvature bounded from below by —K > is given by [TSJ [19], 



K a (x,x;g) < 



V(M) 



+ A{s/2m)- 1 , 



(21) 



where A depends on the diameter of the manifold, d(M), the lower bound of the Ricci curvature 
K, and the volume of the manifold V(M.). Then we similarly obtain 



\U\E)\\ < n 



A 2 r(2) 



r(i/2) 2 J 

+A l / 2 (s(l-u 1 )/m)- 1 / 2 



d s 8 e -^ nm+fi - E y 

- Ul du 2 



1 dui 

KF 2 



v(My/ 2 



V{M) 



T - 2 +A^{s{l-u 2 )/m)-^ 



(22) 



One can even simplify the integrals, that is, the upper bound of the u 2 integral is replaced with 
1 and the square roots of the sums are replaced with the sums of the square roots at the cost 
of getting less sharp bound on the norm of U'(E). Integrating with respect to u±, u 2 and s, we 
have 



\U'(E)\\ < n 



A 2 r(2) 
r(i/2) 2 



i 



V(M) (nm + — E) 2 V(M)V 2 



4A 1 / 2 m V2 7r i/2 r (|)r(i) AmirT( 



1)2 
2 > 



nm + n 



E)2 (nm + [i — E)_ 



.(23) 



In order to get an explicit solution of this inequality, let us put a further natural assumption 
nm + fi — E > u. Then, we find 



\U'(E)\\ < n 



A 2 r(2) 
r(i/2) 2 



K2J V ^ + AmirT(l/2) 2 



1 



V{M)n fiW(M) 1 / 2 
Now if we impose the strict positivity of the principal operator, we obtain 

E gr > nm + /i — n\ 2 F , 

where 



(nm + u — E) 



F 



r(2) 

T(l/2f 



V(M)u 



+ 



4A l/2 m l/ 27r l/2 r( | )r( l : 

uW(M) 1 / 2 



+ AmnT(l/2y 



.(24) 



(25) 



(26) 



Therefore, the lower bounds on the ground state energies for different classes of manifolds (TEJ|) 
and (|25|) are of almost the same form up to a constant factor so the form of the lower bound 
has a general character. From the general form of the lower bounds, we conclude that for each 
sector with a fixed number of bosons, there exists a ground state. However, in two dimensions, 
the ground state energy bound that we have found diverges linearly as the number of bosons 
increases whereas in the three dimensional case it diverges quadratically. Therefore, unlike the 
three dimensional problem, these estimates with our present analysis are good enough to prove 
the existence of the thermodynamic limit in two dimensions. We will now study large n or the 
thermodynamic limit of the model by a kind of mean field approximation, yet this requires a 
more delicate analysis than the three dimensional case. 



4 Mean Field Approximation 



Before applying the mean field approximation, we first clarify one point about this approxima- 
tion, which has not been mentioned in our previous work. It is well known that the residue of 
the resolvent at its isolated pole /i is the projection operator P M to the corresponding eigenspace 
of the Hamiltonian 

1 



p "=-^;f r ,. d£iS(£) - (27 » 

where V ^ is a small contour enclosing the isolated eigenvalue /i in the complex energy plane 
[20J. For the moment, let us consider only the first diagonal element of the resolvent, a(E) and 
choose the contour enclosing the ground state energy E gr which is a well-defined point on the 
real axis thanks to the bound given in the previous section. Then, the above integral of oc(E) 
gives the projection to eigenspace |\E r o)( l Po| corresponding to the minimum eigenvalue. 

From ffl3|) . it is easily seen that the principal operator formally satisfies <f>>(E) = $(£"*). 
As will be shown elsewhere [21], 3>(E) defines a self-adjoint holomorphic family of type A [22J, 
so that we can apply the spectral theorem for the principal operator or inverse of it. Moreover, 
holomorphic functional calculus applies to this family. Most importantly, it can be shown that 
this family has a common dense domain which corresponds to (an operator closure of) the 
domain of H for all sectors. 



E 



\u k (E))(u k (E)\ - doo(E) -L-\ u (E))(u(E)\ , (28) 



where cu(E) (cu k in the discrete case) and \oj(E)) (\oj k (E)) in the discrete case) are the eigenvalues 
and the eigenvectors of the principal operator, respectively. We assume that the principal op- 
erator has discrete, assumed non-degenerate, as well as continuous eigenvalues and the bottom 
of the spectrum corresponds to an eigenvalue. The integrals here are taken over the continuous 
spectrum er($) of the principal operator (for simplicity, we write it formally, it should be written 
more precisely as a Riemann-Stieltjies integral). Due to Feynman-Hellman theorem, we have 



= — ^1 + A 2 J ds s K s (a, a; g) 



e-^ H °-» +E) \u k (E))\ 



-A / ds s 
'o 



e 2 



(Ho-fM+E) 



d 2 x K s (x, a; g)<p g (x)\u k (E)) 



M 



< 



(29) 



by using the positivity of the heat kernel. Note that the operator valued distributions 4> g {x) 
becomes well defined by taking a convolution with the heat kernel in the last term. The bound 
state spectrum corresponds to the solutions of the zero eigenvalues of the principal operator. 
u>k(Eys flow with E due to (1291 . so that the ground state corresponds to the zero of the 
minimum eigenvalue Uq(E) of <&(E). Let us expand the minimum eigenvalue ujo(E) near the 
bound state energy E gr 



cu (E) = u (E gr ) + (E - E gr ) 



du (E) 



dE 



+ . . . = ( E - Egr 



dcuo(E) 



dE 



+ 



(30) 



Using this result and the residue theorem in ( 1271) . we obtain 



H -E gr y l tf{a) 



du {E) 



dE 



\E gr 



\co (E gr ))(u (E gr )\(f>(a)(H Q - E, 



gr I 



(31) 



There is no other pole coming from (H — E)^ 1 near E gr since we assume E gr < nm, and no 
other terms for k ^ contribute to the integral around E gr . Let us assume that the ground 
state eigenvector of the principal operator is 



\ u o(Egr)) 



d 2 x 1 ■ ■ -d?z n _i ip (x 1 , ■ ■ ■ ,x re _i)|xi • ■ -z„_i) • 



(32) 



By using the eigenfunction expansion of the creation and the annihilation operators and their 
commutation relations, we shall shift all creation operators <pL(x) to the leftmost 



Hr, — E 



4>\ (aW ( Xl ) . . . 0t {Xn _ x) = / d 2 . . . d 2 ^ M . . . 1 u yn) 



M r 



X 



/ °° ds e s(H -E + nm) R ^ fl; g) K ^ ^ g) . . . K ^ g) 



(33) 



and all annihilation operators 4> g {x) to the rightmost 

■ ■ ■ Mx^)—^— = [ d 2 g Vl ■ ■ ■ d 2 g y n ds e - S (^-^+nm) 

J^O — hi J M n JO 

x K s (y 1 , a; g)K s (y 2 , a?i; g) ■ ■ ■ K s (y n , x n -i; g)(j) g {yi) ■ ■ ■ <j> g (y n ) , (34) 

which are the generalized versions of the equations we first used in [IT]. Therefore, from the 
equation f[3"Tj) . we read the state vector |\J/ ) 

l*o) = / d 2 g y 1 ---d 2 g y n ^ (y 1 ,...,y n )\y 1 ---y n ) 

= [ d 2 gyi ---d 2 g y n [ d 2 gx l ---d 2 g x n . l - rdse-^ nm - E ^K s ( y(T{lh a;g) 



xK s (y a( 2),Xi;g) ■ ■ ■ K s (y a{n) , x n _u g) V>o(zi, ■ ■ • ,x n -i) 



du {E) 



dE lEar 



-1/2 



\y*{i) ■ ■■y<r(n)) (35) 



where the sum runs over all cyclic permutations a of (123 . . .n). We will now make a mean 
field approximation to this model. In standard quantum field theory, one expects that all the 
bosons have the same wave function u(x) for the limit of large number of bosons n — > oo 
and the wave function of the system has the product form of the one particle wave functions. 
However, due to the singular structure of our problem, the wave function in fl35|) can not have 
a product form in the large n limit. In order to see this, we note that nm — E gr is the crucial 
factor. If nm — E gr = 0(n a ) where a is a positive exponent, we could get a simplification. To 
demonstrate this, let us define nm — E gr — (1 — y[v])x — 2my[v], where y[v] < 1. As we will 
see, x is what we typically estimate, and the variable y[v] is related to the scaled kinetic energy 
functional. Indeed, for compact manifolds we will see that x ~ n1 ^ 2 an d y[v] ~ 0. We scale 
s = s'/[(l — y[v })x — 2my[v}] and as n — > oo, all integrals of the heat kernels are peaked around 
Ua{k)- (This is clear from the property of the heat kernel that K s (x, y; g) — >■ S g (x, y) in the sense 
of distributions as s — > + and also from the stochastic completeness assumption). Then, all 
integrals of x a (i) are 



M 



d 2 g Xi K 8 /[(i_y) x _2my](xi,y a Q + i))lj}o(xi, ...,X h .. . ,X n _i) « lj) (xi, . . . ,y<r(l+i), ■ ■ ■ ,X n _i) 



for I = 1, . . . ,n — 1 as n — y oo and state |\&o) becomes 



(36) 



v&o) « f d 2 gVl ■ ■ ■ d 2 g y n - r ds e- s{nm - E ^K s {y a{1) ,a- g) 

jMn n *e(i-n) Jo 



iPo(ya(2),- ■ ■ ,y<r(n)) \ 7^ \E g 



-1/2 



\y<rQ) ■ ■■y<r{n)) 



(37) 



It is important to note that j^o) is not in the domain of Ho. In order to see this, it is sufficient 
to consider the following term when we calculate (^ol-^ol^o) 

J^d]x J™d Sl e-^ nm - E ^K Sl (x,a;g) J™ ds 2 e -«(™-*W (_-L^ V 2 g K S2 {x,a;g) 



P POO POO / 

/ d]x / d Sl e- sl{nm - E ^K sl (x,a;g) / ds 2 e ~ S2{nm ~ E9r) ( 

J M Jo JO V 



dK S2 (x,a;g) 



(38) 



where we have used the property that the heat kernel satisfies the heat equation — ^V g K s (x, a; g) + 
9Ka ^ a ' 9 ^ = o. After applying the integration by parts to the s 2 integral and using the initial 
condition for the heat kernel K s (x, a; g) — > 5 g (x, a) as s — > + , we find 

/ d]x ^ d Sl e-^ nm - E ^K Sl (x,a;g) 5 g (x,a) 

J M Jo 

poo 

-(nm-Egr) / ds 2 e- S2{nm - E ^K S2 (x,a;g) 
Jo 

POO POO 

= / dsi e~ Sl{nm - Esr) K sl (a,a;g) - (nm - E gr ) / e -^m-E ar ) 
Jo Jo 

/■oo 

x / ds 2 e- s ^ nm - E ^K Sl+S2 (a,a;g) (39) 
Jo 



where we have used the semi-group property of the heat kernel J M d 2 y K Sl (x, y; g)K S2 (y, z\ g) = 
K Sl+S2 (x, z\ g). After the change of variables u = s\ + s 2 and v = s± — S2, we get 

/"OO /"OO 

/ d Sl e siinm-E gr ) K ^ ^ fl . ^ _ ^ _ ^ / &u u e -«(nm-E, r )^^ fl . ^ _ ( 4Q ) 
JO Jo 

The first term is divergent due to the short time asymptotic expansion of the diagonal heat 
kernel for any Riemannian manifold without boundary 



K s (a,a;g) ~ — — V" u k (a, a)(s/2m) k , (41) 

where Uk(a, a, ) are scalar polynomials in curvature tensor of the manifold and its covariant 
derivatives at point a. Similar to the problem with point interactions on manifolds which we 
studied in [13], our problem here can also be considered as a kind of self-adjoint extension since 
the wave function \t>o does not belong to the domain of the free Hamiltonian. The self-adjoint 
extension of the free Hamiltonian extends this domain such that the state is included. 
Although the wave function \l/o is not in the domain of Hq, the eigenfunction ipo corresponding 
to the lowest eigenvalue of $(i?) can be taken in the domain of Hq (see the discussions in |21j). 
As a result, |\l/o) given in ( 137|) is not in the product form in the large n limit, that is, 

„ n 

l*o>^ / d 2 g y 1 ---d 2 gVn ]l* (y k )\y 1 ---y n ). (42) 
Jm " k=i 

The solution takes a kind of convolution of the wave functions in the domain of H with the 
bound state wave function which is outside of this domain. 

Yet, $(-B)'s lowest eigenfunction may be approximated by a product form for large number 
of bosons, that is, 

^o(£i, • • • , x n _i) = u(xi) ■ ■ ■ n(x n _i) (43) 

with the normalization 

|M| 2 = f d 2 g x \u(x)\ 2 = 1 . (44) 
JM 

In our two dimensional problem, for non-compact case, x = O(lnn) as we will see, so the 
product formula ( 1371) may not be a good approximation. Similarly, in three dimensions for the 
noncompact case, x = 0(1). In such a case, w(x)'s are not the wave function of bosons when 
a bound state forms, but related to the correct wave function through (|35|) . In fact, the full 
wave function of the ground state could be read from 

,(n+l)\ \ / h T ,(n+l)\ 



u n) ) r[ 



where we explicitly denote the boson numbers. We call the eigenvector of the principal operator 
in the mean field approximation \u) for consistency of notation with our previous paper |llj . 
In the mean field approximation, the operators are usually approximately replaced by their ex- 
pectation values in this state i.e., (f(x)} ~ /((#)). However, the exact value of the expectation 
value of an operator is given in terms of cummulant expansion theorem if it converges [24] . 
Therefore, we assume that the corrections coming from the higher order cummulants are neg- 
ligibly small and indeed we will see that this assumption is justified for the particular solution 



we will find. Therefore, the expectation value of the principal operator by applying the mean 
field ansatz becomes 



poo 

<j) E [u] = nh [u) - E + fj. + X 2 ds K s (a,a ]g )[e~ s{m -^ -e- s{nho[u]+m - E) ] 

Jo 

POO P 

- nX 2 I ds I d^x d]y K.(z, a; g) K s (y, a; g) u*(x) e --C«*M+*»-s> u ( y ) , (46) 

and 



o Jm 2 



h [u] = I d 2 x ( |Vg 5 X)|2 + m\u(x)\ 2 ) = K[u) + m , 



9 
M 



2m 



where we have taken (n — 1) ~ n for n ^> 1 and K[u] is called the kinetic energy functional. 
Now, we must solve the functional equation (p E [u\ = (giving the bound state spectrum of the 
problem), that is, we solve E as a functional of u(x), and then find the smallest possible value 
of E with the constraint (1441) . One can try to write E as a functional of u(x) from the equation 
4>e[u] = and apply the variational methods to minimize E = E[u]. However, this is a implicit 
function of E and there is no simple way to solve exactly this functional equation since E is a 
complicated functional of u(x). Moreover, we have no explicit expression of the heat kernel on 
any Riemannian manifold to solve E. Nevertheless, it is possible to find a lower bound on the 
ground state energy without applying the variational calculus techniques. 

Since nh [u] and E come together in equation (]4"6"]) . it turns out to be convenient to introduce 
a new variable x = x[ u ] 

X = nh [u]-E. (48) 

Then, the condition (f> E [u] = gives 

-s(m-ju) _ „-s(x+m)l 



POO 

X + V + X 2 ds K 8 (a, a; g) [e~< m ~^ - e 
Jo 

2 

e s( X +2m) _ ( 4Q ) 



f oo 

,2 



n\ z / ds 



d g x K s (x,a;g)u(x) 

M 



Note that the left hand side is an increasing function of x while the right hand side is a 
decreasing function of Xi hence there is a unique solution for %. 

To get a feel for the problem, we consider x as the dependent variable. We now remove the 
X dependence of the right hand side of (149]) by first defining a new dimensionless parameter 
s' = 2m(2m + x)s and scaling the metric cjij = [2m(2m + x)]9ij- Using the scaling property of 
heat kernel in two dimensions 

K s (x, y; g) = a 2 K a 2 s (x, y- a 2 g) , (50) 

and then defining new dimensionless wave function v(x) 

v(x) = [2m{2m + x)}~ 1/2 u{x) , (51) 

all explicit x dependence becomes shifted to the left hand side of (jUJ). The condition 4>e[u] = 
in two dimensions then gives 

(x + A* + X 2 J ds K s (a, a; g) [e^™^ - e - six+m) ] ) 

m 



n\ 2 (2m) 



ds' 



&\x K s r(x,a;g)v(x) 



M 



-s' /2m 



(52) 



Of course, K s >(x,a;g) has now a dependence on \ but let us assume that by varying v(x) we 
can get all possible values of the argument. It is important to notice that the left hand side is 
an increasing function of \ and the right hand side is always positive. Therefore, the left hand 
side is minimum when % = ~~ f 1 - Let us denote the inverse function of the left hand side as 
fi(nU), that is, % = fi( n U[v}). Here U[v] denotes the functional on the right hand side except 
for the factor n. We can express E in terms of the inverse function, 



E = nm + 2mnK[v] + (ni^[v] — l)fi(nll) , 



(53) 



where K[v] = j M d 2 x \V g v(x)\ 2 is the dimensionless kinetic energy functional. Hence unless 
nK[v] < 1, the energy is always bigger than nm + \i. As a result, we see that the interesting 
possibility corresponds to nif [i>] < 1 case. In the analysis for < 1, if we follow the same 

reasoning as in three dimensions, the following integral appears 



ds (1 - e- s{x+2m)/2m )KAa, 



(54) 



However, an upper bound to this integral can not be found by using the diagonal upper bounds 
of the heat kernel given in the previous section because it is divergent for large values of s. 
Therefore, we must develop a different method to handle the two dimensional problem. 

For the case niC[u] < 1, we will again consider the problem for compact and noncompact 
manifolds separately. Using the eigenfunction expansion for the heat kernel [11] and for v(x) = 
J2h^o v (l)fi(, x i d) an d taking the integral with respect to s', we find the right hand side of (152]) 



00 00 
,2 \ ^ \ ^ 



n(2m)A^2Jr 



Z 1= CH 2 =0 



-v*{h)v{h)f h (a;g)f* 2 (a;g) 



(55) 



where fi(x;g) and a\ are the eigenfunctions and the eigenvalues of the scaled Laplacian — V~, 
respectively. Since this is always positive, it is smaller than the following terms by writing the 
zero modes separately 



<n(2m)\'{ |/„(a;s)|> 



E 

(h=0) 



f*(a;g)v*(0)f k (a;g)v(h 



EE 

h?0 h^O 



f^9)v*(k)fi 2 (a;g)v(l 2 



[1 + (T h + (T h 



(56) 



Using / (a;, 



y/V(M®) 



and |f(0)| < 1 we get 



< ra(2m)A 5 



V(M(g)) y/V{M{g)) 



fh(a;g)v(h 



(h=0) 



Oh, 



EE 



fh(a-,9)v*(h)fi 2 (a;g)v(l 2 



[1 + (T h + (T h 



(57) 



I I 



Let us first consider the second term and multiply both numerator and denominator with the 

l-e 

factor cx^ 2 and then apply Cauchy-Schwartz inequality so that we find 

E * ( E hw) 1A ( E JL ) 1A ■ 

where we have chosen < e < 1/2. In order to convert the products (1 + &i l ) 2 d'l~ e in the 
denominator into a summation of them, we use a Feynman parametrization 

1 T(3-e) [\ f 1 6( Ul + u 2 -l)u^ 

dui Ui / du 2 



(i + a h y<t}-* T(2)r(l - e) J J ( Ul (l + a h ) + u 2 a h )^ 

- r(2)r(i - 6) y dni K + ( x Zl )3-,- t 59 ) 

One can express the factor ^ ni+ i p_ e as an integral of s 2 - e e -s ( ui + ''i) due to 

ds s k e~ as , (60) 



a fc +! r(fc + 1 

where 3?(a) > and JJ(fc) > — 1. Therefore, equation ( 159]) becomes 

dm m(l - ui)~ e / — (s72m) 2 - £ e --'(«i+ffi 1 )/2^ . (61) 



(l + oij 2 ^- r(l-e)7 iV 17 7 2m 

Using the eigenfunction expansion of the heat kernel, we have 



^ (1 + ^)^1- r(l-e)7 2m 



x| ^' (a ' fl; ^"v^(I)) 1 ' (62) 

After we transform the above integral to the original metric g and the time s and take the s 
integral by using the diagonal upper bound of the heat kernel for compact manifolds (l2~Tj) . we 
find a upper bound of the above sum 

\f h (a;9)\ 2 ^ AT(2-e) (\ u^ 1 AkT{2 - e) 

^ (1 + ^K" £ " r (!- e ) 'o r(l- e )sinvr e - ^ 

We now come to the crucial point. The upper bound of the first sum in (158jl can be found by 
using the following inequality (also used in |25j). 

< 5 + (^y^ a, , (64) 

where 5 > and < e < 1/2. Note that this inequality applies only to dimensionless variables 
<7i, so that is why we use the scaling transformation at the beginning of the problem as opposed 
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to the three dimensional case. The proof of this inequality is given in Appendix. As a result, 
for the first sum of 0581). we obtain 



2 ^- e < 5 



eN e/l-e 



K\v] 



(65) 



where excluding the zero mode from the sum again gives the kinetic energy functional K[v] 
since a = 0. Then, we finally get 



2n(2m)A 5 



y/V(M(g)) 



£ 

(12=0) 



f h {a-g)v{h 



< 



(1 + *JJ 
2n(2m)A 2 



5 + 



6/1-6 



KM 



1/2 



A7rr(2 



r(l — e) sin7re 



1/2 



.(66) 



Since shi7re > 2e for < ire < tt/2 (a useful inequality: sin 6? > 29 /it for < 9 < it/ 2) and 
nK[v] < 1, the last expression is smaller than 



2 A /n\/2mA 2 



y/(2m + X )V(M(g)) 



6 

n I - 



1 /ev/i-e 



e V5 



" 1/2 


"A7rr(2-e)" 




2r(l-e) 



1 1/2 



(67) 



By choosing arbitrary constants e and 5, we see that the interaction term brings a contribution 
of order 0(n) to the total energy. However, one can even find a better solution to the large 
n behavior of the energy. In order to control the energy as n — > oo, we can assume that the 
parameters e and 5 are sequences in n. Without loss of generality we can assume that e(n) goes 
to zero as n — ?• oo (recall that < e < 1/2). If we want to find a better large n behavior of the 
energy, we must choose the sequence 6(n) such that 



n5{n) 
e 2 {n) 



0(1) 



(68) 



We are now looking for an optimal solution for the energy and tell how fast the sequences e(n) 
and S(n) must change with n. In order to see this, let us write fl67|) in the following way 



2^nV2m\ 2 



y/(2m + X )V(M(g)) 



n6 ( U ) + eT ^yln(6(n)/«5(n))-ln6(n) 

e(n) 



1 1/2 



AttT(2 - e(n)) 
2r(l-e(n)) 



1 1/2 



(69) 



An optimal solution of the sequences can be found in such a way that the exponential term 
goes asymptotically 



: C ^ n ? , ln(e(n)/5(n))— In efn) /^n \ 

e i_ e(n) w ji k ij y ' = 0(mn) . 
This implies that we can choose 



en 



and as a consequence of fl68l) 



<5(n) 



lnra 
1 

rain 2 72 
13 



(70) 



(71) 



(72) 



Therefore, we obtain the upper bound of (I66p for n ^> 1 



2n(2m)A 5 



EJhV_ 
„ (1 



('2=0) 



+ ^1, 



\ 2 \/4mAjre /— — 
< — yn In n . 



^(2m + X )V(M) 



(73) 



As for the last term in (jo7|) . the procedure outlined above is very similar. However, we now 

l-e l-e 

multiply both numerator and denominator of this term with the factor a^ 2 a l2 2 and apply 
Cauchy-Schwartz inequality and get 



EE 



(l + a h +a h ) 



< 



EE^i 



\f h{a;g)\ 2 \fh(g;g)\ 2 



1/2 



(74) 



For simplicity, we can use the following inequality in the second sum 
since cx/ = and obtain an upper bound on ([74]) 



< 



< 



ii^o 



l-e 



§( 1 + ^K 



(l+a tl +*l a ) 2 - (l+a h )(l+a h ) 



(75) 



We again convert the products (1 + criAo] 6 in the denominator into a summation of them by 
using a Feynman parametrization 



1 



T(2 - 



dwi- 



(1 + <7 { X~ e r (! - e ) (1 - ^lY^i + ^) 2 " e 

Then, we rewrite the factor ^ Ui+ ^ p=e using ( 1ST)]) and get 



(76) 



r i 



ds' , 



s72m) 



d«l7 r 

(1 — UiY J 2m 

x J2 e s'(u 1 ^ h )/2m lfh{a .~ g)l 2 



l-e 



(77) 



Using the eigenfunction expansion of the heat kernel above leads to 



du% 



r(l-e)7 (1- Ui yj 2m 



ds' 



(s'/2m) 1 - e e- s ' ui/2m \K s ,{a,a;g) 



V(M(g)) 



(78 



Scaling back to the original variables and using the diagonal upper bound of the heat kernel 
for compact manifolds ( 12T|) . and taking the s integral, we obtain an upper bound of (ITT]) 



A I du 1T ^ 



e-l 



irA 



(1 — u\) e sin Tie 



(79) 



I 4 



The analysis for the first sum in (I75p is exactly the same as before. Therefore, the upper bound 
of the third term in (1571) for n ^> 1 



n(2m)X' 



l + (T h +(Ti 2 ] 



< 2mA\ 2 7reln 



n . 



(80) 



Combining the upper bounds (1751) and fISUl) . we finally obtain an upper of floTj) for n 3> 1 



A 2 



(2m + X )V(M(g)) 



AmireA 



nl/2 



(2m + X )V(M(g)) 



l/2 ln V2 



n 



2mA\ 2 ne 



Inn. (81) 



For the left hand side of ( )52l) . we have 



X + A i <X + / i + -^^/ ds -ftT s (a, a; g) [e" 
'o 



s(m-fi) _ g-s(x+m)j ^ ^Q2^ 



due to the positivity of the heat kernel and x > — A*- Using this result and the upper bound of 
the right hand side of (152]) . we have 



X + V < 



x 2 



(2m + X )V(M(g)) 



n + X 2 — 4m ^^, - NN n 1/2 ln 1/2 n + 2mAA 2 7re Inn. 

(83) 

We will solve \ from this inequality for the large values of n. It is important to notice that 
the right hand side is a monotonically decreasing and the left hand side is a monotonically 
increasing function of x- Therefore, if we find a solution to the above equation, say at x — X*i 
we conclude that x < X*- It is sufficient to find the leading order term of x f° r our purposes so 
X* can be taken as an infinite power series 



1 V2 



(2m + X )V(M(g)) 



X* ~ ain" 1 + a 2 n a2 + 



54) 



with the decreasing power a± > a 2 > ■ ■ ■ of n and ai, a 2 , ■ ■ ■ are the coefficients which can be 
found by substitution. Then we find the leading order term of x 



X< 



A 



\/V{M) 



n 



1/2 



55) 



for n 3> 1. As a result of this, we get the lower bound of the bound state energy E gr for large 
values of n 



E gr ~ nm + jj, 



A 



1/2 



As for the noncompact manifolds, the analog expression of (1551) is 

n(2m)A 2 / / d^h)dfi(l 2 ) * _ v'ihWWhfaflfcfag) 



(86) 



(87) 



where dfi here refers to the formal spectral measure. One can think that the sums are replaced 
with the integrals in the case for noncompact cases and the analysis is basically same as the 
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one for compact manifolds except that we do not have to bother for extracting the zero mode. 
Following the same steps for the analog expression of the compact cases, and using the diagonal 
upper bound of the heat kernel on Cartan-Hadamard manifolds (fT6j) . the ground state energy 
goes like 

E gr ~ nm + /i — 2mCe\ 2 In n . (88) 

Incidentally, these show that K[u] ~ X-^M = y[ v ]x n ~ l ~ y[v\n~^ l ~ aS) where a = 1/2 for com- 
pact and a = for noncompact manifolds. Hence, the cummulant approximation is expected 
to give only the term that we kept, and this shows that the approximation is consistent. 

The solution actually implies an ansatz for the solution which is possible for compact man- 



ifolds. Let us set u 



yJV{M) 

POO 

nm — E + /i + A 2 / ds K s (a, a; g) [ 
Jo 



Then we find for the energy 

-s(m—(j,) s(nm— _B)j 



A 2 



n- 



V(M) \nm - E 



19) 



Now using a scaling for s t—^ s'/ (nm — E) and using the diagonal asymptotic expansion ( I4TI) 
for K s (a,a;g), 



nm 



we find 



E + /x + A 2 



ds' 



2m 
Aits' 



s' (m—fj,) I (nm—E) 



A 2 



n- 



V(M) \nm - E 



m\ 2 , 

nm — E + /i H in 

2n 



nm — E 
m — u 



A 2 



n- 



V(M) \nm - E 



(90) 



(91) 



This has the asymptotic solution as claimed in the compact case. Interestingly, the similar 
ansatz for the three dimensional compact manifolds gives exactly the same type of asymptotics 
for the energy as in two dimensions, which is better than the upper bound obtained of the mean 
field analysis (note that one also has K[u] ~ n -1 / 3 so it suggests a constant wave function as 
an ansatz). 



5 Conclusion 

In this paper, we considered the non-relativistic Lee model on various class of two dimensional 
Riemannian manifolds. It is just a continuation of our previous work [11]. The construction of 
the model and the calculations for lower bound of the ground state energy is almost the same 
as in three dimensions. However, the mean field approximation in two dimensions is not as 
simple as in three dimensions and it required more analysis to predict more precisely the large 
particle limit of the model for compact and non-compact manifolds. 

6 Appendix: The Proof of the Inequality 

In order to prove the following inequality 

x^K S+^j^x, (92) 



where x is a dimensionless variable, let us consider the following polynomial function 



f(x) = 5+ (jj x-x 1 -' . (93) 
We assume that 5 > and < e < 1/2. This function has one extremum point at x* 

*:=(!- e) (-) , (94) 



and this location x* corresponds to the minimum. One can easily see that 

/(x,) = *(l-(l-c) 1 r a ) . (95) 
Since e < 1/2, we have ^ > 1. We can also show that (1 — e) a < 1 — e for a > 1 since 1 — e < 1 

l-e)^<l-e. (96) 



e 

or 



Then, we find 

fM > 5e , (97) 

and it is always positive. Since this is a global minimum point, we obtain f(x) > for all x, 
which completes the proof. 
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